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We investigate integrable fermionic models within the scheme of the graded Quantum Inverse 
Scattering Method, and prove that any symmetry imposed on the solution of the Yang-Baxter 
Equation reflects on the constants of motion of the model; generalizations with respect to known 

■ results are discussed. This theorem is shown to be very effective when combined with the Polynomial 
' ^-matrix Technique (PRT): we apply both of them to the study of the extended Hubbard models, 

for which we find all the subcases enjoying several kinds of (super)symmetries. In particular, we 

■ derive a geometrical construction expressing any gl(2, l)-invariant model as a linear combination 
| of EKS and U-supersymmetric models. Furtherly, we use the PRT to obtain 32 integrable so(4)- 

invariant models. By joint use of the Sutherland's Species technique and 77-pairs construction we 
propose a general method to derive their physical features, and we provide some explicit results. 
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I. INTRODUCTION 



The issue of finding exact solutions for quantum models is strictly connected with the notion of integrability of their 
Hamiltonian H . Integrability actually consists in the possibility of finding a complete set of commuting observables, 
so that the eigenstates of H can be univoquely characterized by the values of the quantum numbers. Within this 
subject, the Quantum Inverse Scattering Method is a very powerful tool providing sets of commuting operators; 
it is based on the Yang-Baxter Equation (YBE) , a functional equation for a C-number matrix (the i?-matrix) , which 
can be thought of as a factorizability condition for the scattering matrix. Such techniques were originally applied to 
study non-linear PDE as well as spin systems, for which some remarkable results have been obtained. As to fermionic 
problems, the traditional approach was to use the Jordan- Wigner transformation to map them into spin systems, 
and then to study the latter by means of the QISM. 

More recently, new methods have been worked out to directly apply the QISM to systems of fermions: in the literature 
two main different approaches can be found to this aim; the first one makes use of algebraic representation theory to 
' find solutions of the YBE, introducing a grading that accounts for the fact that the operators of fermionic systems can 
either commute or anticommute. The second method (fermionic i?-matrix) is quite direct: it consists in dealing with 
t-h ' an operator-valued YBE (instead of a C-number YBE), so that everything is calculated directly in terms of fermionic 
t-H I operators. Despite the remarkable results obtained within each approach, interesting questions are still open, which 
■ we shall investigate in this paper. 

In the first instance, the passage 'C-number' <-» 'fermionic operators' should be clarified in order to use in a comple- 
mentary way the two approaches cited above: in fact the algebraic methods are very powerful in providing solutions 
of the C-number YBE, but their results are not always 'translated' into fermionic operators, so that the Hamiltonian 
and the constants of motion are not easily readable; on the other hand, the fermionic i?-matrix formulation, working 
with operators instead of C-numbers, is more direct but quite cumbersome on a computational point of view. In the 
preliminary sections of this paper we therefore clarify the mutual relationship between these two approaches: sec. II 
and sec. Ill provide systematic methods to represent fermionic operators with matrices and to derive fermion models 
from the C-number YBE. In sec. IV some particular remarks are made on the structure of the graded tensor product, 
which is crucial to this interplay. 

Secondly, a quite important open problem is concerned with the symmetries. Indeed, even if the algebraic methods 
allow to classify the solutions of the YBE according to several kinds of (super)symmetries, no much attention has been 
paid to investigate the effect that these (local) symmetries have on the constants of motion that the QISM determines. 
Although some results have been obtained with the fermionic i?-matrix approach in a particular case Q, just a few 
attempts [|| have been made to find general results on this issue, or to develop them. We point out that knowing the 
symmetries of the constants of motion of exact results is important when developing numerical approaches to study 
perturbations on such exact models; in fact, by fine-tuning the perturbation, one can control what happens to the 
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conserved quantities and realize how many of them get broken. This is crucial because in quantum mechanics it is 
not that easy to have an analogue of the classical Poincare sections for studying the break-up of invariant toruses. 
Moreover, the investigation of symmetries is also very useful to determine the physical features of a model: indeed 
the few exact phase diagrams of interacting fermion models that we know []l3| , have been derived just exploiting the 
symmetry properties. It is therefore quite important, once a given symmetry has been realized to provide interesting 
physical insights, to find all the fermionic models (within a certain general class of Hamiltonians) that fulfill it, in 
order to deduce the properties of further models. Also, it is important to know how many and what kind of free 
parameters the imposition of a certain symmetry allows on a given class of models; in fact, even if the exact solution 
is not possible in general for all the values of the above free parameters, one can select which is the most suitable 
parameter for the investigation of a given physical feature, and see what happens when varing it while keeping the 
others fixed, the symmetry being always preserved. 

We therefore devote all the central sections of our paper to present results on the symmetries. In particular in sec.V, 
concerned with the symmetries of the YBE equation itself, we extend the invariance properties of the YBE to the 
graded similarity transformations, which are the actually useful similarity transformations when dealing with fermionic 
systems. In section VI we prove a theorem assuring that the local symmetries imposed on the i?-matrix directly re- 
flect on the constants of motion; here we shall generalize some known results by considering generic symmetries and 
possible additional signs in the generators of neighboring sites. In sec. VII we investigate the complemantary question, 
namely the conditions under which the symmetries of the Hamiltonian reflect onto the i?-matrix; the Polymonial 
^-matrix Technique (PRT) juj is shown to provide interesting general answers to this issue. Finally, the two results 
of sec. VI and VII are combined together to investigate the specific case of the Extended Hubbard Models (EHM). 
More precisely, in sec. VIII we exploit the matrix representation of sec. II to find all the EHM that fulfill several kinds 
of (super) symmetries (in particular, we provide a geometrical construction allowing to express any gl(2, l)-invariant 
model as a linear combination of the EKS model Jl^Jl5| ] and the U-supersymmetric model ]l7| , p3| | ) ; further ly, in ana- 
lyzing the integrable subcases of each of the above symmetry classes, we realize that for most of them the i?-matrix is 
a first or second order polynomial; hence, we are able to deduce the symmetry properties of the constants of motion 
for such models (e.g. for EKS and U-supersymmetric models). 

Finally, in sec. IX we focus on the so(4) symmetries which are proved to be particularly interesting in condensed 
matter physics. By means of the PRT we find all the models that are derivable by first degree polynomial ^-matrices, 
and through the theorem of section VI we show that the constants of motion are so(4)-invariant as well. All the mod- 
els we find turn out to act as Generalized Permutators; this allows us to exploit the Sutherland's Species technique, 
and to propose a general scheme for finding the spectrum and the ground state phase diagram. The joint use of the 
77-pairs construction (based on the so(4)-symmetry) , leads to determine the eigenstates and the correlation functions 
for these models, even when arbitrary Coulomb repulsion and filling are allowed. 

II. MATRIX REPRESENTATION FOR FERMIONIC OPERATORS 

In the second quantization formulation, the behaviour of a system of fermions on a lattice is described by creation 
and annihilation operators, which are governed by the algebra: 

{Ci,s, c j,s>} = { c i,s, c is>} = kj (!) 

where s is the spin-J label assuming 2 J + 1 possible values (J being a half odd integer for fermions). The local 
d-dimensional vector space Vj associated with the j-th site of the lattice is made up of vectors that are built by acting 
with creation operators on a local vacuum \o)j\ 

Vj = Span ( \a)j = \o), , a = 1 . . . d) (2) 

Here the h^'s are products of m creation operators cj s with different s (m can take values from to at most 

2 J + 1). We shall tipically deal with spin-i fermions, for which a 4-dimensional space is usually involved, so that 

hj = ct.. ; hp — ct 1 ; hj — 1 ; = ct ^ct^.. Due to the anticommutation relations (|l]), the space Vj has an 

intrinsic graduation; in fact Vj — Vj 1 ^ © Vj , where the odd (even) subspace ( Vj) is spanned by those vectors 

that are built with an odd (even) number of creation operators c\ s . Similarly, the space End(Vj) of local linear 
operators on Vj is also graded; odd (even) vectors and operators are also said to have a parity p = 1 (p = 0). 

Moreover, since for any homogeneous 0f\0f ] G End(^) the relation p(Of of ] ) = p{of ] ) + p(of ] ) holds, 
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End(V) is actually a graded local algebra, which we shall denote Aj. Each operator Oj G Aj can be easily given a 
local representation in terms of a d x d matrix by making it act on each basis vector of the space Vj 

of ) \0) j = {O^T p \a) j (3) 

Here (O^)^ are C-numbers, a representing the row and (3 the column of the matrix . It is easily checked that 
this representation is faithful. 

Let us now consider the fermionic problem on the whole lattice; let N be the number of sites. The global algebra 
A is the enveloping algebra of the sum of the local Aj 's. Any global operator can be expressed in terms of linear 
combination of products of local operators: 

O glob € Span ) O ] . . . ojf™ > ) (4) 

where each O?'^ is a generic single-site operator belonging to the local sub-algebra Aj t of the ji— th site, and m can 
run from 1 to TV. In general, the local operators in (Q) need neither appear in the order of the lattice sites nor be 
homogeneous elements of the algebra. 

It is therefore easily understood that the basic tool to set up a matrix representation for operators like is the 
global representation of a single-site operator : since the latter is embedded in the global algebra, the matrix that 
represents it is actually a multi-index (d N x d ) matrix, which can be constructed in principle from the d x d local 
representation determining its action on the local space Vj (see (^|)). However, some caution has to be used in doing 
that, because of the graded structure of both the operators and the vectors. In particular, a graded tensor product 
must be used instead of an ordinary one. 

Therefore, we first of all recall here some notions about multi-index matrices on graded euclidean vector spaces. A 
n-multi-index matrix A is a d n x d n matrix, whose entries are C-numbers denoted by A^'"^° n , each greek index 
running from 1 to d. The upper n-multi-index {a} — (a\, . . . a n ) represents the rows, while the lower n-multi- index 
{(3} — {j3\, . . . p n ) stands for the columns. The couple (a.j,j3j) describes the action of A on the j— th local C d space. 
In each of these euclidean spaces, the canonical basis vectors e a are supposed to be assigned a parity n(a) Ell, so that 
C d is actually a graded vector space. This function ir is in principle completely arbitrary; however we shall choose in 
the following tt = p, where p is the intrinsic parity of the fermionic vectors \a)j of Vj (see(U)). 

The parity of a n-multi-index is defined as p({a}) = p(ot\) + . . . +p(a n ). An n-multi-index C-number matrix A is said 
to be homogeneous with parity p(A) iff for all its non- vanishing entires A^\ one has p({a}) +p{{P}) = const — p{A). 
Given a n-multi-index matrix A and a m-multi- index matrix B (not necessarily homogeneous) , we define their graded 
tensor product as 

= Ag>S«(-ir(W)WW)+ P (W)) (5) 

The graded tensor product is a (n+m)-multi-index matrix. If A and B are homogeneous, then A® s B is homogeneous 
with parity p( A) +p(B). Although the associativity property (A® s B)<3 S C = A <g) s (B ® s C) holds, not any properties 
of the ordinary tensor product can be transferred in general to the graded tensor product; for instance, only if B and 
C are homogeneous one can state that 

{A® S B){C® S D) = {-IY {B ^ C \AC® S BD) (6) 



Let us now come again to the problem of matrix representation of fermionic operators. In order to do that, we remind 
that for fermions a convention has to be specified to define the basis vectors of the global space V g i b = V <S> • ■ • <8> Vn; 
we shall adopt the following one 

\ ai ,a 2 ,...a N ) = f h[ ai) ...h% N) \0) (7) 

where |0) is the global vacuum (defined through c JS |0) = Vi, s), and hj^ is the a^-th of the d operators defining 
the state vectors at the j-th site (see ®)). The action of any global operator O g i b on the global vectors is perfectly 
determined by the algebra ([!]). Its global representation O g i t is defined through: 

O glob |/?i,...,/3jv) = {OgUb) a £;X\«u •••>«*) (8) 
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The key step is to provide a matrix representation for single-site operators; it is straightforwardly realized that, with 
the choice (0) for the basis vectors, the C-number d N x d N matrix representing of as a global operator turns 
out to be: 

0j Q) = I (g> s . . . <g> s I® s O (o) ® s ... ® S I (9) 

j— th 

where is the d x d matrix representing the local action of Oj on the local space (see (^)). Notice that, thanks 
to the associativity property cited above, this definition is unambiguous. 

The matrix representation of an operator O — O^ 31 ^ O^ 32 ' . . . O^ 1 ^ is given by the product of the global representa- 
tions of the single-site operators O^ 3 ^. In particular, if ji — i (i.e. the order in which they appear matches the order 

of the fife's in the definition (0) of the global vectors), then the matrix is simply given by = 0{ ai) s ... <g> s 0^ m) . 
We wish to point out that the matrix representation introduced here also holds for non-homogeneous operators: this 
is very useful because one is not required to decompose a given operator into its even and odd components before 
arriving at its matrix representation. 

A crucial role will be played in the following by Hubbard projectors £j — \a)j j(b\. In the case of a 4-dim. local 
space they are explicitly given by the entries (a— th row and the b— th column) of the following matrix 



/ n, T (l-n.|) 



£s = 



cJ T (l-n it ) 



'in- 1 "H' 

4l c iT n i|(!- n iT) 4lO-~ n iP 

c iT (l - n n ) c n (l - n iT ) (1 - n iT )(l - n n ) 

,t „ „ m J J 



c ji n ji 



" c iT n il 



c il c iT 



c ji n jt 

~ c h n n 
c ji c ji 



(10) 



Each of the above entries is an homogeneous operator with parity p(£j b ) 
enjoy very important properties: 



c b c d 
c 3a ' ° k c 



= 



p(a) +p(b). The Hubbard projectors 

(11) 
(12) 



where [X, Y}± = X Y — (— 1)p( x )p( y ) YX. The matrix representation of the Hubbard projectors will be denoted Ej h ; 
it is given by eqn.(^) where 0^ aj ^ —* E b a , the d x d matrix E b a having vanishing entries except for a 1 at row a and 
column b. The Ej b a 's share the same properties ( |il"| ) and ( p^ ) as the E b h . Using the Hubbard projectors, one can 
express any single-site operator O^ 1 just using its local representing matrix (see ([}])) as follows 



(13) 



III. THE QUANTUM INVERSE SCATTERING METHOD FOR FERMIONIC SYSTEMS 

The Quantum Inverse Scattering Method (QISM) is a powerful tool for studying quantum integrability because it 
provides a set of mutually commuting operators. Within the QISM a key role is played by the so called ^-operator, an 
operator-valued matrix acting on an n-dimensional space termed auxiliary; the /^-operator is local, so that in discrete 
1-dimensional systems one has a matrix Cj for each site j of a chain. The nature of its entries Cj'l determines the 
kind of physical model one is dealing with (i.e. the Cj^ are spin operators for spin systems, fermionic operators for 
fermionic models, etc.); they also depend on two complex parameters which are referred to as spectral parameters u 
and v. When the Cj^ belong to an ordinary algebra Qj, the standard formulation of QISM can be applied as follows: 

if the commutation rules of the ^-operator can be expressed in terms of a n 2 x n 2 C-number matrix R(u, v) through 
the relation R(Cj <8> Cj) = (Cj ® Cj) R (local realization of the Yang-Baxter algebra), then a global realization 
R(T (g) T) = (T ®T)R can be deduced, where Tjy = Cn ■ ■ ■ C\ is the monodromy matrix. From this property a 
set of mutually commuting operators can be obtained Jllfl. The crucial step to pass from the local to the global 
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realization is the relation (£j ® Cj)(Ck®Ci) — (CiCk®CjCi) Vj ^ fc, which stems from the fact that for ordinary 
algebras the local operators referring to different sites commute. 

However, when dealing with fermionic systems, the entries of the ^-operator belong to a superalgebra Q?, CjZ and CkJ 
may not commute, even when involving different sites j and k. Therefore, in order to work out a QISM for fermionic 
systems, a different composition law (other than ®) has to be introduced. To this purpose, it is customary to assume 
that the ^-operator is homogeneous. We recall that in general a matrix A with Grassmann valued entries is said to be 
homogeneous with parity p(A) iff for all its non-vanishing entries p(Afcj) + p({a}) + p({/3}) = const — p(A), where 

p(.<4r^j) is the Grassmann parity of the entry at multi-row {a} and multi-column {/?} (notice that for C-number 

matrices one has p(Ar^j) = and one recovers the definition given in section 11). 

In practice one usually deals with even /^-operators (i.e. p(C) = 0); as a consequence T is even too; in this case the 
required composition law is given by a graded tensor product defined as follows 

^*)c8m#) = A m B \s} (- 1 ) WW)¥W?(W1 (14) 

for arbitrary C/ s -valued matrices A and B (where Q s is the enveloping algebra of the locals Gj 's) . It is possible to 
check that <g> s is associative and that 

{C i ® s Cj){C k ® s C l ) = (dC k ® s CjCi) Vj ^ k (15) 

as it is expected to be for even objects. Thanks to fllq), one can state pi that 



R(u,v) (C(u,w) <g> s C{v,w)) — (jC(v, w) ® s C(u, v)) R(u, v) (16) 

R(u,v) (T(u,w) ® s T(v,w)) = (T(v,w) <g> s T(u,v)) R(u,v) . (17) 

The fact that for fermionic /^-operators a graded tensor product has to be introduced also leads to assuming that the 
R- matrix is even too, i.e. 

P (a)+p(f3)+p( 1 )+p(5)=0 MRf s {u,v)^0 . (18) 



Indeed, this turns out to be very natural when investigating the consistency conditions of eqn.([16|); in fact, one can 
mimic the non-graded case Q and try to derive the Jacobi identities for the structure constants R by reducing a 
multiple product £/ <g> s £/' <g> s Cj'" into its 'reversed' form £/" ® s C/' (g> s £/ through repeated action of eqn.(^6|). A 
compact form is obtained only if R is even too; moreover it can be shown that, despite the fact that graded algebras 
are involved, the obtained relation is non-graded 

(I ® R(u, v)) (R(u, w) (g) I) (I <g> R{v, w)) = (R(v, w) ® I) (I <g> R(u, w)) (R(u, v) ® I) (19) 

This C-number functional equation is known as the ordinary (because no extra signs appear) Yang-Baxter Equation 
(YBE). It is worth stressing that R may depend on the spectral parameters in a completely general way, not only 
through their difference. 

In analogy to what happens for spin models, the above scheme can also be approached backward, i.e. the YBE itself 
can be used to explicitly realize the structure (|l6|), then yielding a set of mutually commuting fermionic operators. A 
peculiar feature of this procedure is that the dimension n of the auxiliary space is taken to be equal to the dimension 
d of the physical local space. Interestingly, one can prove d|J that, starting from an even (see eqn.(|l8|)) solution of 
the YBE ([36]), one can construct an even ^-operator as follows 

C j a {u,v) = {-iy^^Rf s {u,v)E j 5 i (20) 

where the 8j S are the Hubbard projectors. Here a gives the row and (3 the column of the entry. Due to eqn.(|i"l|) the 
entries of two operators Cj and Ck of different sites fulfill 

Cj^(u,v),Ck}(u', v') = Va,/3,7,(S Vu, Vj ^ k (21) 

At the same time the property ( |l2|) allows to show that the Yang-Baxter equation ( |l9| ) is actually equivalent to (|l(|). 
In doing so, one can easily realize that the presence of the additional signs in the definition (En) of the £-operator is 
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crucial; therefore in the graded case, it is also customary to rewrite eqn.(|l^) for the variables R^J = (— l) p ^ p ^ R 1 ^ , 
obtaining the so-called graded Yang-Baxter equation [l|,[To[| 

K%> («» v) R% («, w) R%, (v, w)(-iy (/3')( P («')+Ka")) = (22) 

Starting from eqn.(p^) one can derive a set of commuting fermionic operators. It is worth emphasizing that, since 
R(u, v) is an even matrix, two relations are obtained 

[tr 7jy(u, w), tr TJv(u, w)] = [sir T/v(u, to), strTjq(v, w)] = (23) 

where is the ordinary trace and sir is the supertrace str T = (— 1)p( q )72*. Therefore it might seem that one has two 
kinds of conservation laws; however it is not by developing the equations ( p3| ) powers of the spectral parameters that 
one generates appropriate constants of motion, because neither of them yields local operators; in order to have such 
a property one usually requires to have a couple of values (uq, vq) of the spectral parameters for which the solution of 
eqn.(|l9|) reduces 

r;](u ,v q )=6%S1 (24) 
If this is the case, a straightforward calculation shows that: 

Z(u,v) := (strT A r(u ,wo))^ 1 strT N (u,v) = (trT N (u ,v )y 1 trT N (u,v) (25) 
and that the constants of motion defined as 



^ n = T17 ln - 2: ( u ' w °) 
du 11 



n > 1 . (26) 



are local, in that J n is the sum of operators involving clusters of no more than n + 1 sites Q. The logarithm is taken 
not only to obtain an additive eigenvalue spectrum in performing the Algebraic Bethe Ansatz |2| j2l|j22|| , but also to 
avoid trivially commuting constants. Eqn.(|2^) shows that, although the transfer matrix is usually defined through 
the supertrace 

t(u,v) = strT;si(u,v) (27) 

the ordinary trace plays an equivalent role. However, taking the definition (|2^ ) does simplify the expression of 

the shift-operator, because str TJv(uo, vq) — V^V^ ■ ■ ■ 'Pn-i atj where Vj k is in this case the graded permutator 

Vj k = (~l)P£j^ £k0- Such a simple and compact form is not possible for the trTpj(u, v). A deep argument confirming 

the crucial role of the super-trace will be provided by the investigation of the symmetries (see section VI) . 

Let us consider in detail the first conserved quantity usually interpreted as the Hamiltonian, which is of the form 

H = YliLi with periodic boundary conditions Hn n+i — 'Hni- Each Hj j+i is a two-site Hamiltonian given 

by 

H jj+1 = (-iyW^ +p W d u R$(u,v ) SilSj+ia ( 28 ) 

U — Uq 

Writing down the matrix representation for the Hubbard projectors, it easily verified that the matrix representation 
Hj j+i for the two-site Hamiltonian ( p8| ) is 

H jj+1 =1® S ...® s H^®* ...® s l (H 2sitcs )^ = d u R^(u,v )\ u=u (29) 

3 3+1 

Finally, we also remark that the fermionic operators 1Z and 1Z can be constructed as 

Hjk = (-l) ph){pW+p{S)) R^(u,v)£ J l£4 (30) 
TZj k = (-l)p(f)(pW+P^)R^(u, v) £jl £4 = V a jk Kj k (31) 

Using TZjk and 1Zjk, one can rewrite eqns.([l9|) and ( p2] ) in equivalent forms 
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%23(u,v)'R-l2(u,w)iZ23(v,w) = TZu (v , w) TZ 2 3 (u, w) TZi2 (u, v) 

lZi 2 (u,v)lli 3 (u,w)ll23(v,w) = K23(v,w)K 13 {u,w)1li2(u,v) 



(32) 
(33) 



Also, eqn. (BSf) can be unified into the shorter form 



Hjj+i = d u K"§(u,v ) 



(34) 



IV. SOME REMARKS ON THE GRADED TENSOR PRODUCT 



Throughout the previous sections we used two different kinds of graded tensor product (gtp) ; the former was introduced 
to give matrix representations of fermionic operator (REP-gtp), the latter to adapt the QISM to fermionic systems 
(QISM-gtp). We defined and denoted them differently (<g> s and ® s respectively) to emphasize their different roles. 
Indeed, in order to explicitly convert fermionic operator into matrices or to correctly interpret matrix results in terms 
of fermionic operators, it is worth distinguishing them; since in the literature they are often exchanged and the 
conventions yielding them are not always precised, we wish to comment about that. The graded tensor product <8> s 
(see (||)) fits to represent the fermionic operators when the convention on the definition of the basis vectors is (Q); 

however, had we chosen |ai, a 2 , . . . a/v) d = h^ N) . . . h^\0), we would have got a product like (g> s (see (111)). The 



REP-gtp is therefore strictly related to the convention adopted on the basis vectors. 

On the contrary, the QISM-gtp has different origins: it must be introduced in order to obtain eqn.(|l5|), which could 
not be derived from an ordinary tensor product. The actual explicit form ® s used for the QISM-gtp just stems from 
the additional signs in the definition (|20| ) of the ^-operator; there the indices a and 7 are the rows of R. This choice 
is customary in the literature; however, one could equivalently introduce such extra signs in correspondence with the 
columns, obtaining that the YBE (|l^ ) is equivalent to an expression similar to ([ll]) where £g> s is replaced by ® s . 
Thus, one could in principle use either <£> s or <£> s to define both REP-gtp and QISM-gtp; however, the use of two 
different definitions for the two graded tensor products leads to useful simplifications in the interpretation of the 
matrix results. In fact, the expression ( p8| ) for the two-site Hamiltonian Ti-j j+i contains some additional signs, which 
stem from the definition ( po|) of the ^-operator (so, they are related to the choice of ® s for QISM-gtp); as noticed 
in sec. Ill, when writing down the matrix representation T-ijj+i, the above additional signs cancel out with the signs 
coming from the matrix representation of the Hubbard projectors, provided that <8> s is used as REP-gtp. As a 
consequence, eqn.(p9|) contains no more extra signs. On the contrary, if the same definition ® s was adopted for both 
REP-gtp and QISM-gtp, further additional signs would add to the previous ones, yielding: 

(H 2e ^)7 s = (-1)(p( q )+pM)W)+^» d u R;](u,v Q ) (35) 

U—Uq 

which would contain 'undesirable' extra signs. 



V. THE INVARIANCE PROPERTIES OF THE YBE 



The invariance transformations of a given equation are those transformations that map solutions of the equation 
into other solutions of the same equation. In the case of the YBE the invariance properties are important to un- 
derstand if two models are independently integrable or not. A very important example is supplied by the similarity 
transformations: if a system TL is integrable and we transform the constants of motions through a similarity operator 
Agiob as follows J n — » J' n = Agiob Jn A~ lob , the new operators yield a further integrable system. Nevertheless, it is not 
that obvious in general that the new system can be derived from an i?-matrixQ because local structures may not be 
conserved; however, this is expected to be the case when the transformation is the product of the same local operator, 
i.e. Agiob = A\.. - An, because one is reconducted to study the local change Hjj+i — » (A.jAj+i)Hjj+i(Aj +1 A7 ). 
Actually, for the non-graded case, Kulish and Sklyanin Q proved that this result can still be cast in the language 
of QISM by noticing that the eq. (|l9|) is invariant under similarity transformation of the form A® A; explicitly, if a 



1 for instance, the transformed Hamiltonian may not be of the form Ti! — ^2^L 1 Ti-'a+i 
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n 2 x n 2 matrix R is a solution of ( |l9| ) and 4 is a n x n invertible matrix, then the matrix R' — (A ® A) -R(A ® A) -1 
is also a solution of (|l^). 

However, as to fermionic systems, such a kind of similarity transformations are not suitable; in fact the matrix A still 
represents a local operator Aj, but the matrices representing AjAj+i are built up with the graded tensor product, 
i.e. are of the form A (g> s A; we shall refer to them as graded similarity transformations. The question therefore arises 
whether the YBE equation is actually invariant under such kind of transformations. To this purpose some facts must 
be taken into account: 

i) In section III we pointed out that the Yang-Baxter Equations (|l^) can be set into other forms: we observed that 
when using the variable R — P 9 R ~ where P 9 is the graded permutator {P 9 ) a p J = (-I^WpW Sf - the equation 
( p2| ) is obtained; in non-graded systems a widely used form is achieved rewriting ( |l9| ) in terms of R — PR, where P 
is the ordinary permutator (P 9 )^] = 8$ <fg, arriving at 

Ri 2 (u,v) R 13 (u,w) R 23 (v,w) = R 23 (v,w) R 13 (u,w) Ri 2 {u,v) . (36) 

Although all the above cited forms are equivalent (as long as the matrices are even), it is not trivial in general that 
if the equation in a given form fulfills a certain invariance property, the equation obtained with a change of variable 
fulfills the same invariance property. Now, it turns out that in the non-graded case the invariance property under 
similarity transformation A ® A is actually satisfied from both (|l^) and ([56|); in fact the passages R — > R' — > R' and 
R — > R — > R 1 yield unambiguously the same result, since P{A ® A)P = A® A for any matrix A. On the contrary, 
in the graded case the passage R — > R! — > R 1 has no more the same effect as R — * R —* R', because in general 
P 9 (A®> S A)P 9 ^ A® s A. Thus one has to decide which equation (either (|T9| ) or ([36])) is worth being envisaged. Since 
the Hamiltonian is the derivative of R and not of R (see eqn.(^9|)), eqn.(|19|) seems to be more appropriate. 

ii) the even parity of the R- matrix (which is crucial for constructing integrable fermionic systems from the YBE), 
may not be conserved under such a similarity transformation. Therefore, in order for the transformed matrix R! to 
actually generate a fermionic system, we must take care that p{R') — 0. Interestingly, we shall show that such a 
(physical) requirement is actually the only one that is necessary to prove that the (mathematical) invariance property 
holds. A simple condition to ensure p(R') = is to take a homogeneous A, because in that case A ® s A is always 
even (regardless of the parity of A) and R! is the product of 3 even matrices. Nevertheless, other looser conditions 
are also applicable in principle, depending on both the structure of R and the local A to be considered. 

We can formulate our result as follows: let A be an invertible n x n matrix and R an even solution of (|l9|); if the 
matrix R'(u, v) = (A ® s A) R(u : v)(A ® s A)^ 1 is also even, then R! is still a solution of (|ig|). 

Proof: First of all, it is easily realized that, due to the fact that R is even, the ordinary tensor product in ( |l9| ) can be 
replaced by the graded tensor product <Ss s . Then, the following identities are the key for the proof 

(A ® s A <E> S A)(I <g) s R) = (I ® s R') (A ® s A <g> s A) (37) 
(A ® s A <g) s A)(R <g) s I) = {R! ® s I) (A ® s A <g> s A) (38) 

In order to prove them, one can multiply both equations on the left by (A ® s A® s A) -1 (A <3 S A (® s A), and realize 
that, although for the graded tensor product (A ®) s A ® s A)^ 1 ^ {A^ 1 ® s A' 1 ® s A^ 1 ), it is always possible to write 

(A ® s A ® s A)- 1 = (1 ® s (A ® s A)~ 1 )(A~ 1 ® I ® s I) = (I ® s I ® s A- X ){{A ® s A)- 1 ® s I) 

Moreover, from the property (0), one can obtain that (A® s B)(C® S D) = (AC<E> S BD), provided that B (or C) is even, 
also if the remaining 3 matrices are not homogeneous. Finally, multiplying eqn.(|l9|) by A <®> s A <®> s A on the left, one 
obtains that R! also fulfills eqn.([19|). 



VI. THE SYMMETRIES OF THE CONSTANTS OF MOTION 

In the literature, quite powerful algebraic methods have been developed to classify the solutions of the YBE (|l9| ) 
according to several kinds of symmetries and supersymmetries. An important question is concerned with effects that 
these local symmetries have on the constants of motion J n derived from the QISM (eq.(|26|)). In the literature this 

] and the EKS Q. The standard 
approach for these models is to perform an asymptotic expansion in the (additive) spectral parameter u — v — > oo 
of the monodromy matrix T, and to use the RTT — TTR eq.([l7]) to extract from it the generators of the related 
superalgebra (spl(2, 1) for the t — J and u(2|2) for the EKS). The J n of these particular cases can thus be proved to 
share the symmetries of R. 

In the authors' opinion, however, no sufficiently general approach has been devoted to this subject. To this purpose, 



issue has been examined for some specific fermionic models, such as the t — J |11 
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in this section, by only making use of the existence of the shift-point condition (|24|), we shall prove that the local 
constraints imposed on the R- matrix do reflect onto the whole set of constants of motion. We wish to stress that, 
with respect to the known results concerning specific models cited above, our theorem is more general in that: i) it 
is independent of the dimension of the i?-matrix (i.e. on the dimension of the local vector space Vj); ii) it does not 
make any assumption on the superalgebra to be considered ; iii) it is independent of the specific functional form of 
the i?-matrix with respect to the spectral parameters: in particular, since it does not exploit any expansion on u — v , 
it allows for non-additive i?-matrices, which are recently object of particular interest (see @); iv) it considers also 
the case of 'staggered' operators (the case a = — 1 below). Our results also improve and generalize some arguments 
provided in Q| and with which we shall compare in the following. 

Let us consider a single-site homogeneous operator Xj] as discussed in section II it is possible to write Xj = Xg £j^, 
where (seedl3|)) the matrix X is actually its local representing matrix (||); we obviously have that p(X) = p(Xj). We 
shall show that 

[R(u,v), X® s l + <rl® s X] = a = ±l (39) 

JV 

[J n ,^Xj]=0 (40) 

i=l 

To prove this we shall proceed in four steps: 

Step 1) The condition (39) implies the following relation for the /^-operator entries: 

Cj° X$' + aCj° (-lf™« Xj = Xj £,^(-1)^ + aC& X% Vj (41) 
where the dependence on u and v has been dropped to simplify the notation. 

Proof: It is sufficient to develop ( p9| ) in its matrix entries (e.g. the rows (7, a) and columns (/3,6)) by means of (|^); 
then, by using the fact that Rl" = (—l) p ^ p ^ R^J and that the local dx d matrix X has a definite parity p(X), one 
gets that eqn.(|39|) is equivalent to 

R a p 7 s X$ + oR a p l Xl' (-1)p( x )p W = (_i)pWp(o) XJ, i?;f + aXf, R 5 ^ (42) 

Multiplying this equation by £j S and summing up over 7 and 5, one easily arrives at (^lj) with the help of the identities 

r 01 y — R" 7 V S ' f S anrl Y C a — T) 01 ^' Y a> P S 

3 ft' — - ft /35' ^6 °3 1 Ana ^3 ^3/3 ~ n 06 ^0' C 3 a> ■ 

Step 2) Using ( |4l| ) one can show that a similar relation holds for the monodromy matrix, i.e. 

N N 

T N <$, X% + r N % (-1)p( x MV ]T ° jX 3 = E ^ • ^v|(-l) p(X)p(Q) + <J N T4 Xf, (43) 

3=1 3 = 1 

Proof: This can be done by induction; supposin g th at ( [43| ) holds for chain with N sites, it can be seen that is also 
holds for a chain with N+ 1 sites (for N = 1 eqn. (jf3|) is nothing but eqn. (f4l|) itself). Indeed, it is sufficient to multiply 
eqn. ([43|) by Cn+i^ on the left, sum up over a and use eqn. (ff"l|) with j — ► N + 1; finally, since eqn.(^l|) implies that 

N N 

i=i 3=1 

<*N+l-lNp — [-1) -lN/3 *N+1 



one easily obtains (43) with N — ► N + 1. 
Step 3) Eqn.Q implies that 

if a = +1 



r(u,v) , X 3\ = Vw > u ( 44 ) 

if<r = -l => {t(u,v), EiIi(-l) J '^> =° ( 45 ) 
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Proof: In fact, if p(X) = one takes the supertrace of (f43|), whereas if p(X) = 1 one can take the ordinary trace of 
([43|); in both cases jone obtains equations involving the super-trace of the monodromy matrix, namely ([44|)-(|45|). 
Step 4) Both eqn. (44) and (45) imply that 



N 

~ 1 (u,w)t(v,w) , O -3 Xj 
i=i 



from which eqn. (f40|) is easily deduced due to eqn. (pGj) . 

Proof: For the case a = +1 the proof is trivial; for the case a 

{r-Hw^EiLiC- 1 )^?} = 0, from which 



JV 



= 



-1 one can observe that eqn. 



N 



(46) 



implies that 



[r~V, «0 r(«, to), ^(-1)^-] = 2r- 1 ( U , «;) {r(u, to), ^(-1)^} = 



This concludes the proof. It is worth pointing out that the relations obtained in the cases a = 1 and a = — 1 (see 
eqns.(Q) and ( |45| ) respectively) both involve the supertrace of the monodromy matrix. It would not be possible 
to obtain in general similar relations for the ordinary trace, contrary to what eqn.(p5|) could suggest. This simply 
confirms the role of the supertrace in fermionic systems; nevertheless, a more detailed inspection of the third step 
shows that in the case of even operators Xj, besides eqns.(Q)-(^5|), one also has 



if cr = +1 => 



tr T(u, v) , J2?=i X. 



= 



if <r = -l => {trT(u,v),Y,ti(- 1 ) jx i} = ° v "> 



(47) 
(48) 



meaning that, as long as only even operators are dealt with, the ordinary trace plays a perfectly analogous role as the 
supertrace. 

Finally, it is worth stressing that eqn.(|39]) can be rewritten in terms of the fermionic R- matrix, obtaining the 
equation uZjj+i (u, v), Xj + avfj'+i] = 0. Such form is quite appealing because the symmetries determined by R are 
immediately readable at least for the first constant of motion Jx, i.e. the Hamiltonian, thanks to eqn.([34"|). 

We now want to point out that, for the subcase a = +1, the result ( [40| ) was obtained in 
but from a different kind of local condition on the i?-matrix, which we report here: 



starting not from (J3S|) 



= (-iy( x )iPi a ')+P^))x«,R a '^(u,v) + {-lyW^+P^X^, R*]' (u,v) 



(49) 



where as usual R"g = (— l)" 7 i?^ 
Such an equation was proposed p| just in this form (it cannot be set in a compact form for R) to explicitly suggest 



that, since the ^-operator is constructed with R (see eqn. (20)), the local symmetry constraints should be imposed on 
i?, and not on R. However, it is worth emphasizing that, due to the fact that the constants of motion are derived 
from Z (see eqn.( |25|)) and not directly from the transfer matrix r, the Hamiltonian is directly related to R and 
not to R (see eqn.(|29|)). Indeed the constraint appearing in equation (|3^) is the natural requirement suggested by 
[Hjj+i,Xj + Xj+i] — 0. Anyway, unlike eqn. ([49|) , eqn. (|39|) with a = +1 can by equivalently imposed on R or on 
R. Moreover, it must also be observed that not any operator Xj can be used to implement the condition (|i^) on the 
_R-matrix. In fact, if one evaluates eqn.(|49|) for (u, v) — (uq, vq) and takes into account the fact that for such values 
we have = (— l) ai 5g 81 (see eqn.(pj)), the following relation is obtained 



XjSf 



(l - (-lfW*))) + Xf5} (l - (-1) 



P (x) P (sy 



(50) 



Eqn.(|5Cj|) is identically satisfied when p(X) = 0; on the contrary, when Xj is an odd operator, we can_observe that there 
must exist at least a couple of indices 7, (3 with p((3) ^ p(-f) for which Xj ^ 0; taking a 



S in eqn. (p0[) we obtain that 



in the case cr = — 1 one has to assume that the number of sites is even; this is quite customary in systems that fulfill 
symmetries of the kind ^2iLii~^-Y ( see f° r instance p4). 
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X2 (1 — (— I) 73 ) = 0. This means that eqn.(|49|) only allows for those odd matrices X that have non-vanishing entries 
in odd rows and even columns, which is a somehow 'asymmetric' constraint. This would not allow, for instance, all 
the odd operators of the u(2, 2) superalbegra of the EKS-model (see section VII for more details). 
On the contrary, if one applies the same argument to eqn.(|3^) - or equivalently to eqn.(|42|) no constraints are 
obtained on X, so that any kind of odd operator Xj can be used, as observed above; moreover, one can easily see that, 
if X fulfills the equation (f39|), then X^ does as well, as it is expected to be the case when exploiting the hermiticity 
of the Hamiltonian in [Hj,j+i, Xj + Xj+i] = 0. Therefore the condition ( |39| ) seems to be more general than eqn. ([l9|) . 

Finally, we wish to comment on the physical implication of the above theorem. In doing that we shall anticipate 
some notions on the extended Hubbard models that will be widely treated in the following. The reader non familiar 
with these models can find details in section VII or in the references henceforth given. 

We first of all want to stress that the above theorem (^) => (|l]) generalizes the result obtained by Umeno, Shiroishi 
and Wadati ]?J on the ordinary Hubbard model. Indeed it is known that at half filling the Hamiltonian enjoys the 
so(4) symmetry given by two orthogonal su(2) sectors; the former is given by the spin, and the latter by the operators 
r]~ = J2f = i c ii c ip V + = an d 7 f = J2iLi( n il + n il ~ Using the technique of the fermionic i?-matrix, the 

above authors not only realized that the 7?.-matrix fulfills the so(4)-symmetry, but also showed that all the constants of 
motion enjoy such a symmetry; for the generators rj^ the prove was just based on realizing that {strT^(u, v), n^} = 0. 
In pass, we also remark that, since rj^ are even operators, one could also obtain that {tr Tn(u, v), r^} — 0, as we 
observed above (see eqns(|47|) and (p8[)). A slight modification of the proof supplied in |?J actually confirms that. 
Secondly, it is also worth remarking that, if the i?-matrix is imposed to commute with the number of particles (i.e. 
with Xi = ni), all the J n preserve the total number of particles; then the eigenvectors of the Fock space determined 
by the J n belong to fixed-iV subspaces; this is quite important because it ensures that such eigenstates can also be 
given a first quantization expression in terms of wave functions, which was not obvious a priori since the problem is 
formulated in the non-fixed particle number language of second quantization. 



VII. THE TECHNIQUE OF POLYNOMIAL fl-MATRICES AND ITS ROLE WITH RESPECT TO 

SYMMETRIES 

The Polynomial 7?-matrix Technique (PRT) is somehow complementary (with respect to symmetries) to the theorem 
proved in sec. VI. Indeed the latter exploits the symmetries of the i?-matrix to deduce the symmetries of the the 
constants J n (and in particular of the Hamiltonian); the PRT is a constructive method allowing to look for solutions 
of the YBE starting with a given Hamiltonian of interest; in particular if H belongs to a certain class of symmetry, 
the PRT precises sufficient conditions under which the ^-matrix fulfills the same symmetries as Ti. Using the PRT 
one straightforwardly obtain the result that, if the i?-matrix is a first or a second order polynomial, any symmetry 
imposed on the Hamiltonian immediately reflects onto the i?-matrix. This issue is of great interest for fermionic 
models, in that all the known models (apart from the ordinary Hubbard model, which has very peculiar i?-matrix, in 
that it depends on two spectral parameters) have -R-matrices that are first or second degree polynomials. We shall 
come again to this topic in next section. 

Here we just briefly recall the main aspects of this method, which we shall use in combination with the theorem of 
sec. VI, in order to get further information about the symmetries. We consider for simplicity the case of additive 
R- matrices, i.e. matrices that depend on the spectral parameters through their difference, and search for polynomial 
solutions of the YBE 

R(u) =I+uR^ + ... + ^fl(p) (51) 

pi 

where R^ i — 1, . . . ,p are matrices that do not depend on the spectral parameters; in the case of spin-1/2 fermions 
and 4-dimensional local spaces they are 16 x 16 matrices. 

Inserting the expansion (|5l| ) into the YBE (|l9| ) one obtains a hierarchy of equations for the R^s\ The advantage is 
that such equations are algebraic and not functional. They exhibit some interesting features (see |l4[ ) ; first of all the 
highest degree term (R^ in (jHJ)) must fulfill [^5| the symmetric group equations 

ti$j> f$ R { £ = tig R { £ fl!$ (m) 2 cx I (52) 

In particular this implies that, for first degree polynomial i?-matrices, the Yang-Baxter Equation is equivalent to the 
Symmetric Group equations. 

In addition, the second degree coefficient M 2 ^ is always explicitly given in terms of the first degree one i?' 1 ' as follows 
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£(2) = ^(D) 2 + 7l 



(53) 



where 7 is a C-number. We recall that the first derivative of R(u) with respect to spectral parameter must coincide 
with the 2-site Hamiltonian and therefore R^ must be the representing matrix of the 2-site hamiltonian one is 
interested in. 

From the above discussion, it is easily seen that if R(u) is a first or a second degree polynomial, it always fulfills 
the same symmetries of the Hamiltonian; in fact, if R(u) is of first degree, it is made of the Identity and the 2-site 
Hamiltonian, so that the statement is trivial; moreover, thanks to eqn.(B3T), this property also holds for second degree 
polynomials. 

Combining this observation with the theorem of section VI one can therefore state that every time that an hamiltonian 
7i is reproduced by a first or second order polynomial 7?-matrix, any of its symmetries is shared by the whole set 
of constants of motion J n . In the following we shall apply this sinergic combination to the study of the extended 
Hubbard models; in sec. VIII we deduce the symmetries of the EKS [|l5|,[l2) and the U-supersymmetric model [p"7[~ 
p3j, which are reproduced by a first and a second degree polynomial R respectively. In sec. IX we use the PRT 
to find integrable models that are so(4)-invariant and the theorem to deduce the symmetries of their constants of 
motion. The case of the AAS-model will also be discussed. However, we wish to emphasize the generality of the 
above observations, which can be applied not only to the extended Hubbard models but also for integrable models in 
general. 

VIII. SYMMETRIES IN THE EXTENDED HUBBARD MODELS 

The Hubbard Model is a model of interacting electrons which was introduced to take into account the effect 
of correlations in narrow band insulators; its generalizations - the extended Hubbard models - are envisaged to 
investigate a wide number of physical phenomena such as metal-insulator transitions ]2l| , high-T c superconductivity 
|l5f , quantum wires (2^] , as well as quantum computation Q . We shall consider here the following class of extended 
Hubbard models 

n EHM = - Yl [* - x i n 3,-s + »*,_.) + V-vJ 4,s c k,s + U Y, n iA n u + 

(j,k),8 j 

+ jY, n J n k + Y ^ 4A, S ' c i,s'%s +Y 124,Ai c k,i c k,T + ( 54 ) 

U,k) {j,k),s,s' (j,k) 

O'.fc) O'.fc) J> 

where c j s and Cj s are the usual fermionic operators (see (]l|)). Since electrons are spin 1/2 fermions, the variable 
s can now take two possible values s = {T,!}; the subscript j labels the sites of a lattice A, because the creation 
and annihilation operators are thought of as referring to the set of Wannier functions. Moreover rij jS — cj s c J - s and 
n- = rij* +7i,-i. The symbol (j, k) stands for ordered couples of nearest neighbors in A. In (^) the term t represents 
the band energy of the electrons, while the subsequent terms describe their Coulomb interaction energy in a narrow 
band approximation: U parametrizes the on-site diagonal interaction, V the neighboring site charge interaction, X 
the bond-charge interaction, W the exchange term, and Y the pair-hopping term. Moreover, additional many-body 
coupling terms have been included in agreement with |29(| : X correlates hopping with on-site occupation number, and 
P and Q describe three and four-electron interactions. Finally /j, is the chemical potential. The class (|54|) depends on 
10 parameters and is quite wide; in the literature some even more general class of models have been considered (see for 
instance ]3C| , |3l"| ] ) in which the coupling constants are spin-dependent; however, it can be proved that the Hamiltonian 
([54|) is the most general single-band Hamiltonian that preserves the spin and the total charge, and that is isotropic 
(i.e. 7ij j+i = TLj+ij). It therefore constitutes a fairly general starting point. The su(2)-spin generators explicitly 
read 

S+ = J2 c ]r c Jl S ~ = XX+ c iT S ' = 2^ n it~ n il) ( 55 ) 
3 3 3 

while the w(l)-charge generator is simply given by the operator 
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(56) 



" =2 

3 

In this section we shall investigate for this Hamiltonian several symmetries of the kind envisaged in section VI, i.e. we 
shall impose that J-C BHM commutes with a global generator in the form X = . a 3 Xj, where the single-site operators 
Xj are the generators of a local algebra or superalgebra; the cr = ±1 accounts for possible relative signs between the 
local generators of neighboring sites. 

In this way we find all the mutual relations that the parameters (t, X, X, U, V, W, Y, P, Q, /j,) must satisfy in order for 
these constraints to be fulfilled. To this purpose the matrix representation of fermionic operators developed in sec. II 
turns out to be a very useful tool, since it allows to study the symmetries directly on the C-number matrices represent- 
ing the Hamiltonian and the generators, so that one is reconducted use techniques of ordinary algebra. In particular, 
due to the fact that the % BHM involves nearest- neighbor interaction terms, the global constraint [ H BHM , J2j a ^ ^j] = 
is equivalent to the 2-site form 

[H*™(t,X,X,U,V,W,Y,P,Q,iJ,) , X ® s I + aI® s X] =0 (57) 

where X is the 4x4 local representing matrix of the generator Xj, and H b "" b is the 16 x 16 matrix representing the 
2-site terms of explicitly given in (|5^), where 



hH=hll=fi + V-W h\l = ht\=n + V h\l = hH = h 2 i = hH = (m/2 



hf x = W h% = Y h\l = hjl = -t h\\ = hH = t-2X + X 

h{i=htl=hli = hii = ^ + P+^ + 2V-W hli = hil = fx+^ CM 

h\i = 2/x + 4P + Q + U + 4F - 2W hf A = h\j = -hf A = -h% =t-X 
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(59) 



As observed above, the hamiltonian (^4|) contains 10 free parameters (actually 9, up to an overall multiplica- 
tive i); when imposing a symmetry on H EHM the number of independent parameters obviously reduces; for each 
(super) symmetry that we consider we shall also point out how many of them remain free and comment about it. 

1) 'doubly- occupied sites '-symmetry: It is a su(2) © u(l) © w(l) algebra which preserves, aside the spin su(2) and the 
charge u(l), the following u(l) generator 

K = E(^T-^)K-1-^) ( 6 °) 

3 

In terms of the parameters of ( |5~4| ) this yields just one constraint X = t, leaving therefore 8 free parameters (apart 
from an overall multiplicative factor). This symmetry has been exploited for instance in |32"| , |33"[ to obtain the phase 
diagram of the AAS-model (the subcase X = V = W = Y = P = Q = 0). 

2) so(4)- symmetries: It is formed by two mutually commuting su(2) algebras; six generators are therefore involved: 
three of them are the usual spin components (|5^), while the remaining ones form another su(2) sub- algebra and read 
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tfr) = 12 ° 3 c ^ i c J T V) = 12 ° 3 C J t c i I V z = 2 12 K T + n n - !) ( 61 ) 

j j i 

The case cr = — 1 (which we shall denote by so(4)(_ 1 )) was first considered by Yang S] for the ordinary Hubbard 
model to investigate the symmetry at half-filled band. The case a = +1 (denoted henceforth by so(4)( +1 j) has been 
envisaged for many models |l5| , ^2| , ^3|j3^1 . The generators of both so(4)(_ 1 ) and so(4)( +1 ) are also employed [p4||33| to 
build up states having the property of Off-Diagonal Long Range Order, which is proved to imply superconductivity ]37| ■ 
It can be proved that, in order for (|57|) to hold, the parameters must fulfill the following mutual relations 
so(4)(_ 1 ) (5 free parameters) 

X = 2X V = {W-Y-P)/2 il = Y ~\ Q = -2P (62) 



so(4) /_)_!) (4 free parameters) 

X = t X = V=(W + Y-P)/2 VL= ~ Y ~ l ~2 Q = ~ 2P ( 63 ) 

It must be emphasized that the two fermionic realizations so(4)(_i) and so(4)( +1 ) are physically deeply different: 
for instance so(4)( +1 ) preserves the number of doubly occupied sites (because X = t), while the so(i)( -n does not 
need to. Also, we observe that the number of free parameters that the two realizations allow when ( |57| ) is different 
(respectively 5 and 4); this is because, though it can be found a transformation mapping so(4)( +1 ) into so(4)(_i), 
such transformation does not map the Hamiltonian (|54| ) into itself. 

3) gl(2, 1)- super symmetries: they are formed by 8 generators; the even sector is a su{2) © u(l) subalgebra, made of 
the 3 generators ( |55| ) of the spin and that of the charge (^6|). The odd sector consists of the following four generators 

2 T;(^) =12 cj3 ["t 1 _n J.l) c i,t +l 3n il c j^ 2{. (ct) = (QtiM) 1 

3 

Quw = 12° 3 sJ. ((T) = (q u{ ^ (64) 

3 

In this case the coupling constants appearing in (^4|) become functions of the parameters a and (3 that determine 
such a linear combination. We shall distinguish three cases: 

3a) a 7^ and [3 = : In this case we have 3 free parameters X, U and Q, whereas the other ones must assume the 
values 

X = t W = V = -at Y = -a(t - X) ii = 2at P = 



3b) a = and (3^0: Also in this case the free parameters are X, U and Q, the remaining ones being 

X = t W = -a{t-X) V = -a{t-X) + Q Y = -at fj, = -U P = -Q 

These two cases are related to atypical representations of the superalgebra gl(2, 1) (see 19 2^]). The most interesting 
situation is therefore the following 

3c) and (3 : By denoting b = (3 /a we obtain the relations 

X = t-b(t + a^) X = {l-bf{t + a^) F = | (65) 

W = V= -a[t-b 2 (t + aj)} fi = 2at P = Q = 

where b is a non- vanishing real number. These relations yield two free parameters as long as l + aU/2 ^ 0; in particular 
the subcase characterized by W = V = and 1 + all/ 2 ^ is known in the literature as the U -supersymmetric 
model, and is a 1-parameter model. On the contrary, when 1 + aU/2 — no free parameters are left and the only 
possibility allowed is the so-called EKS-model 

X = t X = U = -2at V = -at W = -at 

Y = -at P = Q = /i = 2at (66) 
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Notice that in this case the values of the parameters are independent of b; this means that the EKS-models commute 
with flsil) for any 6^0 (see the case of u(2, 2) discussed below). 

We also want to remark the relationships between the EKS-models and the U-supersymmetric ones; we remind that 
the U-supersymmetric is the subcase W = V = 0, while the EKS-model has W = —at. First of all it is worth 



rewriting eqns.(65) only in terms of the parameters of the Hamiltonian (i.e. eliminating b). To this purpose the 
couple X-W turns out to be suitable. Since for W ^ at we have two free parameters, whereas for W — —at we just 
have the EKS-models, the allowed values in the X-W plane are those that belong to the two half-planes W/t —a, 
with the addition of the single point (W/t — —a;X/t — 1) representing the EKS model. In the former case eqns.(|65|) 
are rewritten as (we divide by t to eliminate trivial overall factors) 



U/t = 2<7 {- 1 + YTaWjT) X/t= 1 + aW/t P = Q = ° 

(67) 

Y/t = «x (-1 + MA = 2a for: W/t * a 

In the X-W plane the U -supersymmetric model is represented by the axis W — (see fig|l]); the relations for the 
other parameters are the following] 

U/t = 2a ((X/t) 2 - 2X/t) X/t = (X/t) 2 W = V = P = Q = 

Y/t = a ((X/t) 2 - 2X/t) fi/t = 2a (68) 



This pictorial characterization of the U-supersymmetric and EKS model also provides an interesting geometrical 
construction; in fact it can be shown j}8| that any model belonging to the 2-parameter gl(2, l)-class (|6|) can be 
written as a linear combination of the EKS model and an appropriate U-supersymmctric model Tiu BB - In the X-W 
variables this property can be expressed as follows: 

U gl{2>l) (X, W) = [t + aW] H Us3 - aW Ueks (69) 

The interesting fact is that the appropriate model Ti.u ss can be determined in a geometrical way; given any gl(2, 1)- 
invariant model (characterized by a certain point P of coordinates (X, W)) it is sufficient to trace the line from P to 
the point E — (X/t = V.W/t = —a) (which is the EKS-model), and determine the point A at which the above line 
intersects the axis W = (see hg.|l|). This point precisely represents the TLu ss -model we need in eqn. (|69|) . Its coupling 
constants can be determined from its ascissa Xa by means of the relations (|68|). 

Eqn.(|69|) is quite important in that it explicitly shows that the study of gl(2, l)-invariant models with W ~ —at 
actually consists in treating the U-supersymmetric model tiu„ as a perturbation of the EKS-model; on the contrary, 
in the limit W w 0, it is the EKS-model that acts as a perturbation on the known model Tiu BB ■ 



4) so(5) -symmetries: it is an algebra made of 10 generators; we shall consider the realization in terms of fermionic 
operators introduced in |3q] 



A 23 


= 25+ 


A 32 


= 2S- 


A 22 = 


2S Z 


A 14 


= Vd 


A 41 




A u = 


2 V * 


A 12 


= E^t 

3 


A 21 


j 


A 13 - 


-E 

3 



3 A 



a J c 



31 



A 31 



3 



a J c 



(70) 



This is a subcase of the algebra so(4)( +1 ) defined above; the obtained relations among the parameters only allow 
the EKS-models (|66|). The latter are therefore the only models of the extended Hubbard class (54) that enjoy the 
so(5)-symmetry; actually, there are other 1-d single-band models commuting with (|70| ) (see |35f| ), but they do not 



3 The first relation of (g8|) shows that, in spite of the name, the U-supersymmetric model is more suitably parametrized by X 
than by U. Indeed U is a continuous function of X (a parabola), while X is not even a single- valued function of U : this drawback 
both forces to consider only a sub-part of the parabola, and may also lead to non appropriate choices of the parametrization 
(like in [fill) that wrongly let think of discontinuities at U = which actually do not exist. 
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belong to the class (|5J) since they are not isotropic. Even if such models exhibit a superconducting behaviour, it 
should be observed that a more realistic description of the superconducting transition is reached when assuming that 
the Hamiltonian is isotropic and the wave function breaks the lattice symmetries (otherwise there is no break-up in 
fact). 

5) u(2, 2) -super symmetry: It is made of 16 generators; the even sector contains the 3 generators of the spin and the 
3 ones of the so(4)( +1 ) algebra, as well as the identity and the doubly-occupied sites generator (|60|). The odd sector 
contains the following eight generators 

3 3 

QU°)=J2 ajn 3l C 3l QU°)=J2 ajn 3l C n 2t;(a) Cj;(«r) ( 71 ) 

3 3 

The commutations with its generators only allow the solution ( |66| ) of the EKS models. If one considered the algebra 
so(4)(_]/) instead of so(4)( +1 ) in the even sector, no extended Hubbard model j54| ) would be found. 



We wish to make some remarks about the EKS-modcls, defined by eqn.(|6q). The original EKS-model introduced 
in |1| corresponds to the case a = +1 of cqn.(|6^); however, we have unified a = ±1 in the definition since they 
often share the same kind of symmetry. Here we precise the conditions for this to ha 
shown that the 2-site term of the EKS-modcls arc respectively a graded permutator (15 



jen. First of all, it is easily 
(the case a = +1, denoted 

V 9 ), and an oppositely-graded permutator (the case a = —1, denoted P 9 ): their 16 x 16 matrix representations are 
{ p %8 = {-l) p{a)pW SfS? and (P 9 )^ = (~1)(p( q )+ 1 )(p(' 3 )+ 1 ) SfS? respectively. It can also be proved that, for any 
homogeneous operator Xj, 

r%X,Vl X k ^.V,7^ ( 1 v -V- . (72) 

Using ( [72|) it is easy to realize that, if X is the local representing matrix of Xj 

[P 9 , X® S I + I® S X] = VI even and odd (73) 



whereas 

[Pa , X® S 1 + 1® S X] = VXeven (74) 
[Pa , X ® S I-I<Z) S X] = VXodd (75) 

These relations imply that, when investigating the commutation of the Hamiltonian with generators of the form 
J2i <J " 'Xj, both the EKS-models are found to fulfill the commutations as long as additional a 3 appears in odd gener- 
ators (see for instance the cases of gl(2, l), so(5) and u(2,2) symmetries); in particular V 9 is obtained if a = +1, 
while V 9 when a = — 1. Indeed odd generators with a = +1 can be mapped into those with a = — 1 through the 
transformation Cj — > (— 1) J Cj (which leaves unaltered any even generator). 

On the contrary, when the sign a 3 appears in even generators, we always obtain V 9 for the subcase a = +1, but not 
V 9 in the subcase a = — 1 (see for instance the cases of the so(4) symmetries). Indeed the case of symmetries with 
respect to even generators with a = — 1 is quite peculiar; in particular the study of the generators rf^s ( |6l| ) with 
a = — 1 is physically intriguing for reasons of energetic stability of the ground state [^9| ; we shall investigate more 
deeply this aspect in section VIII. 



Here we want to briefly recall which are, among the several classes of symmetry found above, the subcases that 
(up to now) have been proved to be integrable. As discussed in the previous sections, the integrability of a nearest- 
neighbor Hamiltonian TL is proved whenever it is possible to find a matrix R(u, v) which solves the YBE ( |l9| ) and 
such that its first derivative with respect to the spectral parameter equals the representation of the local hamiltonian 
Hii+i (see eqn. (§§)). 

Focussing on the symmetries and the supersymmetries considered above, we remind that as to the 'doubly occupied 
sites'-symmetry, the i?-matrix has been found for 96 models (subcases of which are known models). Within the 2- 
parameter class of gl(2, 1) invariant models the i?-matrix has been found for the subcases of the EKS-models jlj] and 
the 1-parameter U-supersymmetric models |l9|-pl]|. For the so(5)-symmetry and the u(2, 2)-supersymmetry the two 
EKS models are the only extended Hubbard models ( |54| ) that fulfill them, as observed before. So, their ^-matrices 
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are the ones cited just above. As to the so(4)-symmetry, the ^-matrix has been found for the half filled Hubbard 
model @]. 

As anticipated in sec. VIII, there is an intriguing feature shared by the ^-matrices of all the integrable models cited 
just above (apart from the Hubbard model, which has a non-additive R- matrix): they are all polynomials in the 
difference of the spectral parameters; more precisely, the EKS-models have first degree polynomial _R-matrices, while 
the U-supersymmetric have a second degree one. Even if they are not usually written in the literature as polynomials 
(because the form in which they are given depends on the method used to find them as solutions of the Yang-Baxter 
equation), it is easy to realize that such i?-matrices can be cast into polynomials through a mere multiplication by 
a scalar function tp or redefinition of the spectral parameters (see [Q))- This fact means that all the known models 
(apart from ordinary Hubbard) can be found by means of the Polynomial i?-matrix Technique discussed in sec. VIII. 
In particular we therefore have that all the constants of motion of the U-supersymmetric model enjoy the gl(2, 1)- 
supersymmetry. As to the EKS models (|6^) the conclusion is still more general; for the case a = +1 we have the graded 
permutator, and thus, using eqn.(f73|) we deduce that the constants of motions J n of this models fulfill any global 
symmetry J2j %j\ f° r the case ct = — 1 we have an oppositely graded permutator and therefore (due to eqns.(|74|)-(|75|)) 
we have that the constants of motion J n commute with Ylj %j f° r anv X even, and with ^2j(—l) J Xj for any X odd. 
Furtherly, first degree polynomial R- matrices for 96 models have been found in [ p5| ; all these Hamiltonians are shown 
to fulfill the 'doubly occupied symmetry' su(2) © u(l) © u(l); now we can deduce that all the constants of motion of 
these models fulfill such symmetry as well. 

Among the symmetries envisaged in this section, the case left to be considered on the point of view of integrability 
is the one of so(4)-invariant models. This class is very large (contains a number of free parameter), and in spite of 
the fact that it is quite interesting on a physical point of view, it has been investigated only partially up to now. We 
shall therefore devote the whole section IX to it. The PRT will be applied together with the theorem of section VI to 
find so(4)-invariant models and to deduce the symmetries of their constants J n . 



IX. INTEGRABLE MODELS WITH SO(4)-SYMMETRIES 

The class of so(4)-invariant extended Hubbard models (see eqns.(|6^) and (|63|)) is relevant in physics because, by 
exploiting its properties of symmetry, some interesting features on the eigenstates of these models can be easily de- 
duced, such as the form of the correlation functions, which exhibit a Long Range Order. It is then important to have 
as many exact results as possible for such class. We recall that it is a rather wide class, since so(4)(_i) allows 5 free 
parameters and so(4)( +1 ) 4 free parameters. 

Among the known so(4)-invariant models, the most famous is perhaps the ordinary Hubbard model at half filling 
p4j (which corresponds to the subcase X = W = Y = P = Q of so(4)(_ 1 j, see eqs.(|62[)). The ^-matrix for this 
model was found in Q; also it has been proved that its constants of motion fulfill such symmetry as well. Later, 
the EKS-models and the AAS model were also found, for which the phase diagrams were derived; they enjoy the 
so(4)( +1 )-symmetry. Moreover, in the model recently studied by Alcaraz and Bariev within the Coordinate Bcthe 
Ansatz, an so(4)(_ 1 )-invariant subcase (r/ = and e = +1 in pl[) can be identified; it corresponds to the 1-parameter 
subclass X = (1 - sm6)t ; Y = -W = t cos 9 and P = of eqs.(|§). 

Motivated by the physical interest in the so(4)-symmetries, in this section we shall find further models enjoying such 
symmetry. More precisely, we shall make use of the polynomial i?-matrix technique to find all the integrable so(4)- 
invariant extended Hubbard models that are derived from a first degree polynomial R(u) in the spectral parameter 
u. We consider here both so(4)( +1 ) and so(4)(_i). The number of such models turns out to be 32, 16 with symmetry 
so(4) ( -_ 1 ) and 16 with so(4)( +1 ); the EKS and AAS models are shown to be among them. In providing the ^-matrix 
for all of them, we shall guarantee their integrability; moreover, by exploiting the result of the theorem proved in 
section VI, we will be able to show that not only the Hamiltonian but the whole set of constants of motion of these 
models are so(4)-invariant. All these models have the natural properties of preserving the total magnetization and 
charge and of being isotropic, since they belong to the class ( p^ ) by construction. 

In order to show how to find them, we have to use as R^ the matrix (|5^) of the extended Hubbard models, with 
the additional prescriptions (^) and ( |63| ) of the so(4)-symmetries, and impose that R^ fulfills eqns.(p2|). We give in 
the following the values of their coupling constants, as well as the i?-matrices that yield them. In order to simplify 
the notation we divide the 32 models into four subgroups according to vanishing or non-vanishing pair hopping and 
exchange amplitudes; they are denoted Ti.( a \ TL^ h \ H> c > and HS d ^ respectively. The result are presented in a compact 
form: the case a — +1 gives the models that are invariant under so(4)/ + u whereas a = —1 corresponds to the 
so(4)(_ 1 )-symmetric cases. 

I s * subgroup: H (a) = 8 models with Y ^ andW ^0 
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X = t X = (l-a)t U = 2 Sl t V = Sl t W = -s 2 
Y = crs 1 t P = -(s 1 + s 2 )t Q = 2(s 1 +s 2 )t n = -2s!t 



(76) 



)//(/ 



subgroup: TiS b ^ = 8 models with Y ^ and W — 



) rd 



X = t X = (1 - a)t U = 2s x t V = (si + s 2 )t W = 

r = crsit P = -(s 1 +2s 2 )t Q = 2(si + 2s 2 )t n = -2stt (77) 

subgroup: Ti.^ — 8 models with Y = and W =/= 

X = t X = (l-a)t U = Asxt V = sit W = s 2 t 

Y = P = (-2si + s 2 )f Q = (4si - 2s 2 )t f i = -2 Sl t (78) 

4 th subgroup: = 8 models with Y = and W = 

X = f X = (1 - a)t U = Asxt V = (si + s 2 )* W = 

F = P=-2(si + s 2 )t Q = 4( Sl +s 2 )i ^ = -2*i* (79) 

For the Hamiltonians the subcases with s 2 = — Si and a = +1 are the two EKS models, while the subcase 
s 2 = — Si = 1 and a = — 1 is the model proposed in fjlf with 77 = = and e = +1. No U-supersymmetric model 
appears, since one always has X ^ X 2 . The case HS d ^ with s 2 = —s\ and a = +1 is the AAS-model [ B^ , g3| . 
The models ©-(0)-©-© 

can be obtained from the following i?-matrices 

R(u) =l + u (if 3 ( £ tM + Sl tl) (80) 

where i? 2 _sites with i = a,b,c,d is the 2-site matrix (^9|) of coupling constants (^6|) , ( |77|) , (fr^) and ( |79|) respectively 
^. Since the above ^-matrices are first degree polynomials in u, we have by construction that [R(u) , X so ^ ® s I + 
<7 1 ® s X so ( 4 ) ] = where X so ^ is the local representing matrix of the single-site generator = cj^ct-p i.e. the 4x4 
matrix (see sec. II). Thus, thanks to the theorem of section VI, we can deduce that 

j 

We wish now to discuss how to deduce some physical properties of the models found above. We shall provide here a 
general scheme of techniques that allow to obtain information about the spectrum, the phase diagram of the ground 
state and the behaviour of the correlation functions for the 32 models presented above. Details on the application of 
the method to some specific case are given in separate publications (4^^^]. At the end of this section we provide as 
an example some results obtained for two models of group 2. 

In the first instance we shall describe the method to derive the spectrum of the above models; such method is based 



on the observation that the Hamiltonian of the above 32 models (|7q) -j- (79) can be proved to be of the form: 



WW = ^11,.,.; (81) 

3 

(up to the trivial additive term s\ tl) where IL,j + i is a 2-site generalized permutator. A 2-site generalized permutator 
(GP) acts on two neighboring sites exchanging some couples of states but leaving unchanged some other couples; more 
explicitly: 

\a) 3 \(3) j+1 = 9° ap |/3), \a) j+1 + 9 d ap [a^ |/3>,+i (82) 



4 As it is clear from eq.@, the A-matrix (gg) yields H (a) , H (b) , H {c) and H (d) up to an additional energy shift sitl, which 
obviously changes nothing to the integrability of such models. 
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where 0° ^ and 6^ are two discrete- valued (0 , + 1 or - 1 ) and 'complementary' (| 0°^ | = 1 — 1 6 d p | ) functions that identify 
the exchange / non-exchange respectively; notice that in both cases an additional sign is possible (if 9°,8 d = — 1); 
moreover, 0°^ = 9^ a , due to hermiticity. These functions 9 d p and 9° a ^ completely determine the specific kind of GP 
one deals with. 

In pass, it can be proved p5| that within the class of the extended Hubbard models ( |5 4 |) there are 96 (and only 96) 
models that can be cast in the form (|8l]), with a suitable GP; the 32 models ([76]) -f- (\l9j) that we have presented here 
constitute the subclass of all the so(4)-invariant GP, and it is also possible to see that they cannot be mapped into 
other GP through any unitary transformation. 

Remarkably, recognizing that a Hamiltonian has the structure ( |8l| ) allows to exploit the Sutherland's Species Technique 
jl0|,[ll]], which actually consists in regarding a GP as an ordinary permutator between the so-called Sutherland's species 
(SS), the latter being groups of the d local physical states. The groups of states are not determined by the original 
Fock Space but uniquely by the structure of the Hamiltonian; the SS therefore relate to the dynamical processes 
involved in the model. In the case of the extended Hubbard models (154), the physical states per site are 4, and 
thus the SS are always no more than 4. Following Sutherland's notation |41|, each species can be classified as either 
'fermionic' (i* 1 ) or 'bosonic'(-B) (odd and even in j!o|, respectively), according to the sign of the non-exchange part 
(9 d = +1 for B, while 9 d = —1 for F). In particular, within the 32 so(4)-models, one can distinguish the following 
types of SS: 



I s * subgroup F A 

2 nd subgroup F 3 

3 rd subgroup F 3 

4 th subgroup B 2 



B 2 F 2 ; B A 

BF 2 ; B 2 F ; B 3 (83) 
B F 2 ; B 2 F ; B 3 
BF ■ F 2 



where B n F rn characterizes a model with n bosonic and m fermionic SS. Different models that are recognized to be of 
the same SS-types also share the same structure of Coordinate or Algebraic Bethe Ansatz equations; one is therefore 
reconducted to known problems, as the possible different signs of off diagonal terms (6 d ) are not expected to alter 
the structure of such equations. For the cases F p or BF P the Bethe Ansatz equations have been investigated |4l]] 
within the framework of the Coordinate Bethe Ansatz. The cases F p and B p have been examined in |42| within the 
QISM. Some others cases [BF 2 and B 2 F 2 ) have the same algebraic structure as known models (t — J [[ll 16| and 



EKS J12| ) |15j respectively). For the BF type the spectrum, up to constant terms, is that of a spinlcss fermions on a 
chain: 

e({n l }) = -2j2co8(^)n l (84) 

1=0 

where m = 0, 1 are quantum numbers and L the length of the periodic chain. 

The whole spectrum can be found in all these cases. However, it must be realized that the actual degeneracy of each 
eigenvalue depends on the way the Sutherland's species are realized in terms of physical species in each specific model. 
This in turn enters the calculation of the partition function, determining different physical features (for the discussion 
in case of BF models see [Q ) . 

Within the spectrum, the ground state is particularly worth of interest, because these Hamiltonians are expected to 
well describe materials that exhibit peculiar physics at low temperatures. For the F p case, it can be shown |DJ that 
the minimum of the energy, whose eigenvalue reads 

e = l--j o d X — , (85) 

is reached at equal densities of all fermionic species. Interestingly such result, when applied to the F 3 and F A models 
in (^), implies that the ground state in these cases always contains doubly occupied sites, which are expected to 
model short coherence length pairs. 

Also, again as far as the ground state is concerned, it has been shown (i^] that Sutherland's theorem (originally 
formulated f4l| for ordinary permutators) can be extended to the generalized permutators, and it is thus possible to 
assert that, in the thermodynamic limit, the ground state energy (per site) eo of a B n F m problem is equal to that of 
a BF m problem. The latter types become therefore the only relevant ones as to the issue of determining e . 

We also want to stress another important general feature of the 32 models ([76]) ~ ((79|) ; their Hamiltonians H.^\ aside 
the so(4)-symmetry, also preserve the number of doubly occupied sites (see section VII), since they all have X = t. 
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Hence, by adding to each 7iW (i = a , b, c, d) further terms of on-site Coulomb repulsion and chemical potential with 
arbitrary coupling constants, 



n >(i) =n (i) +u Y / n jr n jl +»Y,( n n+ n jO > ^ = a,b,c,d , (86) 

3 3 

the obtained models H' ^ still preserve the number of spin-up, spin-down electrons separately, as well as the number 
of doubly occupied sites. Such property helps to diagonalize the Hamiltonian Tt' W within each subspace of given 
eigenvalues of S z , if and IC; this is interesting from a physical point of view, because one can investigate how the 
features of the model change when tuning either the parameter U or the filling n (which can be expressed in terms of 
fj,); the former is somehow an intrinsic energy unit for these models, and is actually the most meaningful parameter for 
the systems where strong electronic correlations are involved, as pointed out by Hubbard himself in his original paper 
p6[ ; the latter is also relevant in that it can model the degree of hole doping in the material. The two parameters U 
and n are expected to drive the transitions between physically different phases (Metallic, Insulator, Superconducting). 
The so(4)-symmetry of the part TL^ of eq.dSrjl) implies that ]^] 

[W' (i) ,< ) ] = (f/-f)< ) ■ (87) 

Eq. ([§71) is crucial in order to implement the so called 77-pairing mechanism for the construction of eigenstates of TC W : 
indeed, once a reference eigenstate \ref) otH'^ is given, further eigenstates ^(o-)) can be found by applying rjf ff s on 

I re/), namely ^(o-)) = ijf{a)) m \ re f)- This property allows to reconduct the calculation of the correlation functions for 
the eigenstates ^(o-)} to those of the reference states \ref), which are in general easier to be computed. Within the 
class of the above 32 models, meaningful reference states can be found: the vacuum |0) is always one (the eigenstates 
1^(0-)) constructed on it are called pure 77-pair states); the eigenstates \U = 00} of the U — 00 Hubbard model |3(| are 
also reference states whenever W — V = 0; similarly, for the cases in which W — V — ±i the eigenstates \t — J) of the 
t - J-modcl @0 are reference states. In these cases the eigenstates ^(o-)) constructed on \U = 00) or on \t — J) 
are referred to as mixed 77-pair states. The correlation function for these reference states have been investigated in 



the literature (see for instance |lq,p3 34 1) 



The techniques outlined above are quite general and apply to all the 32 so(4)-models <$7(\) -7- (|79|) . As already ob- 
served, all the models Tt^ of a given B n F m -type share the same spectrum equation structure. However, this does 
not necessarily imply that the related W'W have all the same physical features; indeed the presence of the on-site 
Coulomb repulsion U has two main effects: i) it can change the eigenvalues and their degeneracy; ii) it can lead, 
according to how the SS are actually realized in terms of physical species, to quite different shapes of the ground state 
phase diagram as a function of U, even for models of the same type (see for instance ^o) for the -B-F-type). 

We just wish now to provide here a concrete example of the use of the scheme proposed in this section, by presenting 
results on the ground state of the models ( |86] ) with = TL^ (subgroup 2), and s\ = — S2 = — 1; the coupling 
constants of these two models explicitly read X — t; X = (1 — <r)t ; Y = at ; W = V = ; Q = —2P = 2t; U and /z can 
be taken as arbitrary, as observed above. The part W 1 ' is a CP for 3 Sutherland's species, namely F = { f), | j)}, 
Bi = |0), and B 2 = \ IT}- Using the extension of Sutherland's theorem one can derive the structure of the ground 
state phase diagram as a function of U and n; details of this derivation can be found in |40| . The result is presented 
in fig.|. 

The structure is the same for both models a = ±1. Four regions can be recognized: in region I the ground state 
is made of doubly occupied and empty sites; the eigenvalue of the energy (per site) reads: eo = (U/2 — l)n; the 
eigenvectors are pure ry-pair states; the behavior of the two-particle reduced density matrix 

_ ni{y(*)\4Ai c 3i c n 

^^'■J = Ttii Kft \ 

m(*(o-)|*(<T))m 

in this region is the following |15|,|33| 



5 To simplify the notation, we omitted any subscript a in denoting Ti^ and Tl'^h we precise that, in eq.(^), one has to take 
( res P- V(Li)) when H (t) in TC ' <l) contains a = +1 (a = -1); see eqns.@-r(||). 
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(pa)ij " e^ 1 ^ ^ n n (l - n n ) (89) 



where nf± is density of doubly occupied sites. As eq. (|88j) shows, in region I one has a Off-Diagonal Long- Range-Order 
in this state, which implies that the model is superconducting in this region. It is also possible to see that, with respect 
to the AAS model, where the pure 77-pair phase is degenerate in the momentum of pairs |33f| , here the pair hopping 
term selects a specific momentum. This momentum (which is in principle observable through neutronic spectroscopy) 
is strictly related to the kind of so(4)-symmetry allowed by the model; in particular, for the case of so(4)f + -n (i.e. 
a = +1 in eqns.([77|)) the selected momentum is q = 0, while for the case of so(4)(_ij (i.e. a = —1 in eqns.(fnj)) the 
pairs have momentum q = n. The EKS model has 0-momentum pairs (see Jl5|). 

In region II all the possible types of sites (empty, singly and doubly occupied) in the ground state are present; the 
energy reads 



^ = (l-|)(iarccos(i4)-n)-^l-i(l-| )2 , (90) 

and the eigenstate is mixed 77-pair, constructed on the \ref) = \U = 00) states. The correlation function behaves in 
this case Q as: 

"n ^X)"'* ((i - - (91) 

where n s is the density of singly occupied sites. This region is again superconducting and is particularly interesting 
in that it survives up to positive values of U (which are expected to be more physically meaningful, because the 
electronic on-site interaction should be repulsive, even if partially screened by the phononic effective attraction). This 
feature is shared by the other known models like the AAS and the EKS ones; however it must be observed that, 
around half-filling, the region II of our models raises up to U max — 6t, the highest value of TD exactly known models. 
Finally, region Ill-a is metallic; the eigenvalue is eo = — 2ir~ 1 sin(7rn), while the eigenstates are those of the \U = 00) 
model; region Ill-b is the particle-hole transformed of region Ill-a (eo = +2tt~ 1 sin(7m)). At half-filling it is possible 
to show that a charge gap A = U — 6t exists, which makes the model an Insulator. 

This is the picture of the ground state; for excited states one cannot make use of Sutherland's theorem and is re- 
conducted to the spectrum of a whole B 2 F case. On the physical point of view, when the temperature is turned 
on, thermal fluctuations are expected to break the Long-Range-Order of the superconducting regions, according to 
Mermin- Wagner theorem. Nevertheless a Quasi-Long-Range Order can survive, because the decay of (p2)i,j can pre- 
serve a tail over macroscopically observable distances. 

Finally, we want to briefly comment on the mutual relationship between the 32 so(4)-invariant models US 1 '. Within 
each of the 4 subgroups, there are transformations mapping some models into others: in particular a model THS^ 



characterized by (s\, S2, cr) in its parameters - see eqns.([76|)-^(79)- is connected to the one with (— si, —S2, cr) through 
the transformation Cj iS — > (— ) l Ci yS ; analogously, a model withfsi, S2, cr) in its parameters can be mapped into that 
with (si, S2, — cr) through Ci jS — > [1 — (1 — (— Y)ni i - s ]ci :S . Also, the transformation c^j — > (— l) 8 c] , maps the subgroup 
H [h) into the H (c) . 

An exhaustive classification of all these transformations is out of the purposes of the present paper, and will be 
treated in a forthcoming paper. Here we just want to emphasize the following aspects. In the first instance the above 
transformations act differently on neighboring sites; this implies that they are not graded similarity transformations, 
because (in matrix representation) they are of the form A® s B, and not A® s A (see section V). Indeed, by exploiting 
the matrix representation developed in section II it is possible to see that none of the above models is connected 
to any other through a graded similarity transformation; therefore, as far as the i?-matrix is concerned, they are all 
independently integrable. Secondly, it is also possible to show that the models with different values of the parameter 
Q cannot be connected by any kind of transformation (not only similarity). Thirdly, it should be pointed out that, 
even if some transformation maps a given model into another, it is not that obvious in general that the ground state 
of the former is mapped into the ground state of the latter. This confirms that, also for what concerns their physical 
stability properties, all the above models deserve a deep interest. 

X. CONCLUSIONS 

In the preliminary part of this paper we have clarified the relationship between the two approaches to fermionic 
integrable systems, by providing a systematic method to pass from fermionic operators to matrix representation 
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and from the C-number YBE to fermion models. Further, we have proved general results about the symmetries in 
integrable systems: we have treated the symmetries of the YBE, showing that it is invariant under graded similarity 
transformations. Also, we have proved, under quite general hypothesis, that the symmetries imposed on the i?-matrix 
directly reflect onto the whole set of the constants of motion J n that the QISM determines. The complementary 
question whether the symmetries of the Hamiltonian reflect on the ii-matrix has been widely developed with the 
PRT. These two general results are shown to be very effective when combined together; in particular we applied both 
of them to the study of the extended Hubbard models (EHM). By making use of the matrix representation discussed 
above, we found all the EHM that fulfill different kinds of symmetries and supersymmetries ('doubly occupied sites', 
so(4), gl(2, 1), so(5), u(2, 2)), also including in the generators possible relative signs between two neighboring sites. 
In particular, for the 2-parameter subclass with gl(2, 1) supersymmetry, we have provided a geometrical construction 
expressing any gl(2, 1) model as a linear combination of the EKS-model and a U-supersymmctric model. Furtherly, 
by showing that most of the integrable known model (such as the EKS and the U-supersymmetric and others) are 
reproduced by first/second degree polynomial 7?-matrices, we deduce the symmetries of their constants of motion. 
Finally, focusing on the case of so(4)-symmetries (whose interest in condensed matter has been discussed), we exploited 
the PRT to find all the EHM that are derivable from first degree polynomial ^-matrices. The symmetries of the 
constants of motion of these models have also been discussed. Further, by making use of the Sutherland's species 
technique, we have proposed a general scheme to derive the spectrum and other physical properties of these 32 
models; in particular it has been observed that, by means of Sutherland's theorem, the ground state energy can be 
determined for all of them. Finally, thanks to the so(4) symmetry, the addition of arbitrary Coulomb repulsion and 
chemical potential terms to these integrable models can be used to implement the r/-pairs construction. We proposed 
an example of concrete application of the described techniques to two of the 32 models, for the ground state of which 
we have explicitly given the eigenvalue, the eigenvector and the pair correlation functions in the different region of 
the phase diagram (fig.|j). 

We also point out that all the extended Hubbard models investigated above can also be used as 'bulk systems' to 
which one can add appropriate boundary interaction terms; this would lead to new results in the context of models 
with Kondo impurities (see for instance H|). 
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FIG. 1. The X-W plane to describe the 2-parameter class of gl(2, l)-invariant models. The coupling X stands for correlated 
hopping amplitude, while W is the exchange term. The allowed values are the two half planes W/t > —a and W/t < —a, 
with the addition of the single point E = (1, — a) that represents the EKS models (a = ±1). The 1-parameter subclass of 
U-supersymmetric models are represented by the X-axis (W = 0). 

Any gl(2, l)-invariant model (point P) can be written as a linear combination of the EKS model and a particular 
U-supersymmetric model, whose parameters can be geometrically determined by tracing the line from P to E, and finding the 
ascissa Xa of the point A in which it intersects the X-axis. 
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FIG. 2. Ground state phase diagram of the two models X = 1; X = (1 — a); Y = —a; P = — 1; Q = 2 (subgroup 2, eq.(^)) 
from ref.[41]. The model exhibits an insulator-superconductor transition at n = 1, for £/ c = 6. The dashed line is the EKS 
model, and the dotted line corresponds to the AAS model. 
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